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CHAPTER  I 


INTRODUCTION 


The  U.S.  Army  Corps  of  Engineers  is  responsible  for  the 
operation  of  hundreds  of  multiple  purpose  reservoirs  in  addition 
to  maintenance  of  hundreds  of  miles  of  non-reservoir  projects 
(e.g.,  levees  and  navigation  channels).  Management  of  reservoir 
releases  for  water  quality  can  be  analyzed  to  determine  the 
operation  with  any  one  of  the  numerously  available  reservoir 
computer  programs  [WRE  1969a,  HEC  1972,  U.S.  Army  1977,  HEC  1978, 
Loftis  1980].  With  river  water  quality  programs,  the  impact  of 
specified  reservoir  releases  can  be  evaluated  at  downstream 
points  of  interest  [HEC  1978]. 

The  problem  with  using  single  project  models  is  the 
difficulty  of  coordinating  releases  among  projects  which  impact 
on  a  single  location.  This  is  particularly  obvious  in  Figure  1 
where  the  operation  of  both  reservoirs  A  and  B  impact  on  the 
amount  and  quality  of  water  at  City  A  (i.e.,  control  point  3). 

As  the  system  is  expanded  further  downstream,  the  computations 
necessary  to  provide  a  best  operation  of  reservoirs  A  through  D 
for  control  point  7  obviously  require  a  comprehensive  system 
approach. 

The  "HEC- 5 ,  Simulation  of  Flood  Control  and  Conservation 
Systems,  Appendix  on  Water  Quality  Analysis"  computer  model  [HEC 
1986]  has  been  developed  specifically  for  evaluating  the  type  of 
problem  shown  in  Figure  1.  The  model  is  capable  of  evaluating  a 
reservoir  system  of  up  to  ten  reservoirs  and  up  to  thirty  control 
points.  The  model  will  determine  a  best  system  operation  for 
water  quantity  and  quality;  evaluating  operational  concerns  like 
flood  control,  hydropower,  water  supply,  and  irrigation 
diversions.  Changing  needs  and  natural  inputs  can  also  be 
accommodated . 

The  HEC-5  water  quality  routines  can  be  used  to  determine 
the  quality  constituents  available  with  the  best  system  water 
quantity  operation  or  alternately  the  best  water  quality 
operation.  The  model  can  analyze  water  temperature,  up  to  three 
conservative  and  three  non-conservative  constituents,  dissolved 
oxygen  and  phytoplankton.  Optional  computation  intervals  from 
hourly  to  monthly  are  available.  Graphical  post-processor 
capability  can  be  interfaced  through  other  available  software. 

This  training  document  provides  guidance  on  the  application 
of  the  HEC-5  computer  program  to  a  typical  water  quality  study. 
The  purpose  of  this  training  document  is  to  familiarize  the  first 
time  user  of  HEC-5  with  the  procedure  to  follow  for  collecting, 
assembly,  and  manipulating  water  quality  input  data.  The 
optional  types  of  executions  and  the  proper  interpretation  of 
results  are  also  discussed  at  some  length.  The  author  conveys 
many  significant  items  not  normally  discussed  in  a  users  manual 
or  even  in  short  course  lectures.  These  items  resulted  from 
experience  gained  by  completing  several  studies  with  this  water 
quality  model. 


Following  the  procedures  in  this  document  will  help  the 
reader  apply  HEC-5  to  routinely  encountered  problems  involving 
evaluation  of  water  quality  conditions  in  existing  and/or 
proposed  multipurpose  reservoir  systems. 

The  HEC-5  water  quality  model  is  new  and  therefore  a 
research  tool  until  it  has  been  successfully  applied  to  numerous 
practical  problems.  The  HEC  would  appreciate  your  comments  and 
observations  which  could  be  added  to  this  report  or  to  the  users 
manual  regarding  experiences  with  the  application  of  HEC-5.  As 
desirable  improvements  are  identified,  modifications  will  be 
made . 


Additional  assistance  in  understanding  HEC-5  is  available  by 
contacting  Mr.  Willey  at  (916)  551-1748  or  (FTS)  460-1748. 


Figure  1 

TYPICAL  RESERVOIR  SYSTEM  SCHEMATIC 
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CHAPTER  II 


MATHEMATICAL  MODEL 

The  mathematical  model,  "HEC-5,  Simulation  of  Flood  Control 
and  Conservation  Systems,  Appendix  on  Water  Quality  Analysis," 
has  been  used  to  analyze  the  Sacramento  Valley  [Willey  1985], 
Kanawha  River  Basin  [Willey  1986],  and  the  Monongahela  River 
Basin  [Willey  1987]  Reservoir  Systems.  The  computer  program 
Users  Manual  [HEC  1986],  and  several  technical  papers  [Duke  1984, 
Willey  1982,  1983,  1984,  and  1987]  adequately  document  the 
details  of  the  model  concepts  and  the  input  description,  so  only 
a  brief  overview  is  provided  in  this  chapter. 

The  HEC-5  water  quality  computer  program  is  composed  of  a 
flow  simulation  module  (HEC-5A)  and  the  water  quality  simulation 
module  (HEC-5Q) .  They  are  an  integrated  package  with  feedback 
capability  between  the  two  modules.  Each  module,  the  gate 
selection  routine  and  the  flow  alteration  option  are  described 
below. 

Flow  Simulation  Module 

The  flow  simulation  module  was  developed  to  assist  in 
planning  studies  for  evaluating  proposed  reservoirs  in  a  system 
and  to  assist  in  sizing  the  flood  control  and  conservation 
storage  requirements  for  each  project  recommended  for  the  system. 
The  program  can  be  used  in  studies  made  immediately  after  the 
occurrence  of  a  flood  to  show  the  effects  of  existing  and/or 
proposed  reservoirs  on  flows  and  damages  in  the  system.  The 
program  should  also  be  useful  in  selecting  the  proper  reservoir 
releases  throughout  the  system  during  flood  emergencies  in  order 
to  minimize  flooding  as  much  as  possible  and  yet  empty  the  system 
as  quickly  as  possible  while  maintaining  a  balance  of  flood 
control  storage  ("balanced  pool")  among  the  reservoirs. 

The  above  purposes  are  accomplished  by  simulating  the 
operation  of  a  system  of  reservoirs  of  any  configuration  for 
short  interval  historical  floods  or  for  long  duration  nonflood 
periods  or  for  combinations  of  the  two.  Specifically  the  program 
may  be  used  to  determine: 

a.  Flood  control  and  conservation  storage  requirements  for 
each  reservoir  in  the  system. 

b.  The  influence  of  a  system  of  reservoirs  on  the  spatial 
and  temporal  distribution  of  runoff  in  a  basin. 

c.  The  evaluation  of  operational  criteria  for  both  flood 
control  and  conservation  (including  hydropower)  for  a 
system  of  reservoirs. 

d.  The  expected  annual  flood  damages,  system  costs,  and 
system  net  benefits  for  flood  damage  reduction. 

e.  The  system  of  existing  and  proposed  reservoirs  or  other 
alternatives  that  result  in  the  maximum  net  flood 
control  benefits  for  the  system  by  making  simulation 
runs  for  selected  alternative  systems. 
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Water  Quality  Simulation  Module 

The  water  quality  simulation  module  was  developed  to 
simulate  temperature,  as  well  as  three  user-selected  conservative 
and  three  user-selected  non-conservative  constituents.  The  model 
also  allows  dissolved  oxygen  to  be  simulated  if  the  user  selects 
either  carbonaceous  or  nitrogeneous  oxygen  demanding 
constituents.  An  option  for  phytoplankton  evaluation  is  also 
available . 

The  water  quality  simulation  module  accepts  system  flows 
generated  by  the  flow  simulation  module  and  computes  the 
distribution  of  all  the  water  quality  constituents  in  each  of  the 
reservoirs  in  the  system  and  their  associated  downstream  river 
reaches.  The  reservoirs  may  be  in  any  arbitrary  parallel  and 
tandem  configuration. 

The  water  quality  simulation  module  also  selects  the  gate 
openings  for  reservoir  selective-withdrawal  structures  to  meet 
user-specified  water  quality  objectives  at  downstream  control 
points.  If  the  objectives  cannot  be  satisfied  with  the 
previously  computed  "balanced  pool"  flows,  the  model  will  coirpute 
a  modified  flow  distribution  necessary  to  better  satisfy  all 
down-stream  objectives.  With  these  capabilities,  the  planner  may 
evaluate  the  effects  of  proposed  reservoir-stream  system 
modifications  on  water  quality  and  determine  how  a  reservoir 
intake  structure  should  be  operated  to  achieve  desired  water 
quality  objectives  within  the  system. 

Each  reservoir  is  assumed  to  be  a  control  point,  in  keeping 
with  the  concepts  used  in  the  development  of  the  flow  simulation 
module.  Additional  control  points  may  be  placed  in  the  stream 
system  below  the  reservoirs  at  stream  confluences  and  any  other 
desirable  locations. 

Computational  time  steps  from  hourly  to  monthly  are 
optional.  The  model  is  limited  to  simulations  of  one  calendar 
year. 

The  reservoirs  are  represented  by  a  series  of  one¬ 
dimensional  horizontal  elements  such  as  those  shown  in  Figure  2. 
Each  horizontal  element  is  characterized  by  an  area,  thickness 
and  volume.  In  the  aggregate,  the  assemblage  of  layered  volume 
elements  is  a  geometric  representation  of  the  prototype 
reservoir.  This  one-dimensional  representation  has  been  shown  to 
represent  adequately  water  quality  conditions  in  many  deep,  well 
stratified  reservoirs  by  Eiker  [US  Army  1977],  Baca  [1977]  and 
Water  Resources  Engineers  [1968,  1969a,  1969b]. 

Each  horizontal  layer  is  assumed  to  be  completely  mixed  with 
all  isopleths  parallel  to  the  water  surface  both  laterally  and 
longitudinally.  External  inflows  and  withdrawals  occur  as 
sources  or  sinks  within  each  layer  and  are  instantaneously 
dispersed  and  homogeneously  mixed  throughout  each  element  from 
the  headwaters  of  the  impoundment  to  the  dam.  It  is  not 
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Figure  2 

GEOMETRIC  REPRESENTATION  OF  A  STRATIFIED  RESERVOIR 
AND  MASS  TRANSPORT  MECHANISM 


possible,  therefore,  to  model  longitudinal  variations  in  water 
quality  constituents.  Simulation  results  are  most  representative 
of  conditions  in  the  mam  reservoir  body. 

Vertical  aavection  is  governed  by  the  location  of  inflow  to, 
and  outflow  from,  the  reservoir.  Thus  the  computation  of  the 
zones  of  distribution  and  withdrawal  for  inflows  and  outflows  are 
of  considerable  significance  in  operation  of  the  model.  The  WES 
withdrawal  method  [Bohan  1973’  is  used  for  determining  the 
allocation  of  outflow.  me  Lu-vier  -.'.low  allocation  method 
[Debler  1959]  is  used  for  the  placement  of  inflows. 

Vertical  advection  (physical  movement  of  mass  due  to 
continuity  balance)  is  the  net  interelcment  flow  and  is  one  of 
two  transport  mechanisms  used  in  the  module  to  transport  water 
quality  constituents  between  elements.  Effective  diffusion  is 
the  other  transport  mechanism.  The  effective  diffusion  is 
composed  of  molecular  and  turbulent  diffusion  and  convective 
(physical  movement  of  wa*  »r  due  to  dor.sit}  !  *'*  tnfoi  lityj  nixing. 

Wind  and  flow-induced  turbulent  dir  fusion  and  convective 
mixing  are  the  dominant  components  of  effective  diffusion  in  the 
epilimnion  of  most  reservoirs.  Tn  quiescent ,  we  1 1 -strati f ied 
reservoirs,  molecular  di' lusiou  may  be  a  significant  component  in 
the  metalimnion  and  hypolimmon.  For  deep,  well-stratified 
reservoirs  with  significant  inflows  to  or  withdrawals  from  the 
hypolimnion,  flew  induced  turbulence  in  the  hypolimnion 
dominates.  For  weakly  stratified  reservoirs,  wind  induced  or 
wind  and  flow  induced  turbu’ent  diffusion  will  be  the  dominant 
component  of  the  effective  diffusion  throughout  the  reservoir. 

One  of  two  methods  may  be  selected  by  the  user  to  calculate 
effective  diffusion  coefficients.  For  shallow  weakly  stratified 
reservoirs,  the  wind  controlled  mixing  [HEC  1978]  method  is 
appropriate,  while-  the  stability  method  j  HEC  1978]  is  more 
appropriate  for  deeper  well  stratified  reservoirs.  Both  of  these 
methods  have  been  shown  in  numerous  applications  to  adequately 
represent  the  mixing  phenomena  for  heat  and  dissolved  water 
quality  constituents  when  nroperly  applied. 

The  stream  system  is  represented  conceptually  as  a  linear 
network  of  segments  or  volume  elements  as  shown  in  Figure  3. 

Each  element  is  characterized  by  length,  width,  cross-sectional 
area,  hydraulic  radius,  energy  slope,  Manning's  n,  and  a  flow  and 
depth  relationship.  Flow  rates  at  stream  control  points  are 
calculated  within  the  flow  simulation  module  using  any  one  of  the 
several  programmed  hydrologic  routing  methods.  within  the  flow 
simulation  module,  incremental  local  flows  (i.e.,  inflow  between 
adjacent  control  points)  are  assumed  to  be  located  at  the  nearest 
control  point. 

Within  the  water  quality  simulation  module,  the  incremental 
local  flow  may  be  divided  into  components  and  placed  at  different 
locations  within  the  stream  reach  (i.e.,  that  portion  of  the 
stream  bounded  by  the  two  control  points).  A  riow  balance  is 
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GEOMETRIC  REPRESENTATION  OF  STREAM  SYSTEM  AND 
MASS  TRANSPORT  MECHANISM 


used  to  determine  the  t ..  ow  it  < men*'  boundaries.  Any  flow 

imbalance  (i.e.,  the  d  i  f ference  i  n  th^  flow  at:  the  upstream 
control  point  plus  all  tributary  inflows  and  the  flow  at  the 
downstream  control  point)  is  distributed  uniformly  to  the  flows 
at  each  element  boundary.  Once  interelencnt  flows  are 
established,  the  depth,  surface  width,  a . ; ;!  -toss  sectional  area 
are  computed  at  each  clement  boundary  either  by  using  the  input 
flow-depth  relationship  or  by  .ssum no  r.  rr.  1 1  depth. 

Gate  Selection  Rout  ir.  - 

Once  the  desired  resejvoir  .  e’en  ;■>  and  *• h<-*  target  water 
quality  to  meet  down:  t  mam  ne.  las  been  t  o muted ,  the  cate 
selection  algorithm  uetarminev  *h : t.  h  ports  should  be  open  and 
what  flow  rate  should  pass  through  each  open  port  in  order  to 
maximize  a  particular  function  of  the  downstream  water  quality 
target  concentrations.  Solution  of  this  problem  is  accomplished 
by  using  mathematical  optimization  techniques.  The  objective 
function  is  r  •  • :  ••  •  •  •  •  ■  m..-m  -  •  :■  virM  •  n.ri  1  i  t es 

subject  to  Vd  ■  x.o.:  ;  t:  ;  •  - ■  ir.-.i  :  vidua  1  ports. 

The  reservoir  intake  bt  *-uc:ture  con  have  up  to  two  wet  wells, 
containing  up  to  ei  qh*  •  ■".=  .  :n,  -nd  3  'T'-oi  control  outlet. 

It  is  assumed  that  rc  .  e  ■$,  t  n.  ;>-gh  :,y  r  t  these  ports  (including 

the  flood  control  outlet)  leave  the  nservoir  through  a  common 
pipe.  At  any  given  time,  only  one  port  in  either  wet  well  and 
the  flood  control  outlet  may  be  operated.  Hence,  the  algorithm 
provides  flows  through  three  r  "  t s  at  "c  T . 


The  HEC-5  model  also  provides  for  releases  through  an 
uncontrolled  spillway.  These  releases  are  not  a  part  of  the  gate 
sexection  algorithm,  but  the  wafer,  quality  of  the  spillway 
releases  are  considered  by  the  ,  v  e  select  ion  algorithm. 

The  algorithm  proceeds  ny  considering  a  sequence  of 
problems,  each  representing  a  different  combination  of  open 
ports.  For  each  combination,  the  optimal  allocation  of  total, 
flow  to  ports  is  determined.  The  combination  of  open  ports  with 
the  highest  water  quality  index  defines  the  optimal  operation 
strategy  for  the  time  period  under  corn ideration . 

There  are  four  different  types  of  combinations  of  open 
ports.  For  one-port  problem,  a’  of  the  flow  is  taken  from  a 
single  port  and  the  w-iter  qua  I  i  ..y  ..icy  is  computed.  For 
two-port  problems,  ;  omb \ cat  :mi  i .  of  one  port  in  each  wet  well  and 
combinations  of  each  port  with  the  flood  gate  are  considered. 

For  three-port  problems,  combinations  of  one  port  in  each  wet 
well  and  the  floodgate  are  considered.  The  total  flow  to  be 
released  downstream  is  specified  externally  to  the  gate  selection 
routine,  but  if  the  flow  alteration  option  is  selected,  then  the 
flow  can  be  treated  as  an  additional  decision  variable  and  the 
flow  for  which  the  water  quality  index  is  maximized  is  also 
determined . 


For  each  combination  of  open  ports,  a  sequence  of  flow 
allocation  strategies  is  generated  using  a  gradient  method,  a 
gradient  project icn  method,  or  a  Newton  projection  method  as 
appropriate.  The  value  of  any  flow  allocation  strategy  is 
determined  by  evaluation  of  a  water  quality  index  subject  to  the 
hydraulic  constraints  of  the  system.  The  sequence  converges  to 
the  optimal  allocation  strategy  for  the  particular  combination  of 
open  ports.  These  problems  are  solved  very  efficiently  by  using 
mathematical  optimization  techniques  that  take  advantage  of  the 
problem  structure,  namely  a  quadratic  objective  function  with 
linear  constraints. 

Flow  Alteration  Routine 

The  flow  alteration  routine  is  designed  to  change  the 
reservoir  releases,  computed  by  the  flow  simulation  module,  to 
better  satisfy  the  stream  control  point  water  quality  objectives. 
Timing  of  intervening  tributary  inflows  are  considered.  Second 
order  effects,  such  as  reaeration  and  external  heating  due  to 
increased  or  decreased  stream  surface  area,  are  not  included. 

The  calculation  procedure  for  the  flow  alteration  option  is 
as  follows: 

1.  The  relative  mass  of  a  water  quality  parameter  being 
simulated  that  needs  to  be  added  to  the  flow  at  the 
control  point  (for  those  constituents  below  the  target) 
or  reduced  in  the  flow  at  the  control  point  (for  those 
constituents  above  the  target)  is  computed. 

2.  The  average  reservoir  release  concentration  is  computed 
for  all  reservoirs  for  which  the  constituent 
concentration  in  the  releases  is  greater  than  the 
target  concentration  at  the  control  point  of  interest 
(for  those  constituents  below  the  target)  or  for  which 
the  constituent  concentration  in  the  releases  is  less 
than  the  target  at  the  control  point  of  interest  (for 
those  constituents  above  the  target) . 

3.  The  total  dilution  flow  requirement  is  then  computed  by 
dividing  the  result  of  step  1  by  the  result  of  step  2 
to  provide  the  total  flow  release  needed  to  bring  the 
constituent  concentration  at  the  control  point  of 
interest  to  the  target. 

4.  The  flow  is  then  apportioned  to  the  reservoirs  capable 
of  bringing  the  control  point  constituent  concentration 
to  the  target  in  proportion  to  the  flows  originally 
computed  for  those  reservoirs  by  the  flow  simulation 
module. 


Thus  the  flow  alteration  requirement  can  be  computed  for 
each  control  point  and  for  each  constituent.  The  various 
computed  flow  rates  are  then  combined  by  using  the  coefficients 
of  a  linear  programming  objective  function  and  the  deviation  of 
the  respective  constituent  concentrations  from  the  target 


concentrations  at  each  respective  control  point. 

Once  the  flow  augmentation  requirement  is  determined,  the 
flow  simulation  module  is  recalled  and  the  computations  for  flow 
and  water  quality  are  repeated  for  the  final  results. 

Summary 

HEC-5  model  is  capable  of  simulating  the  water  quality 
effects  of  the  operation  of  a  system  of  reservoirs.  Each 
reservoir  may  be  operated  to  satisfy  a  number  of  objectives, 
including  flood  control,  low-flow,  hydropower  production,  water 
supply  and  water  quality  control.  The  water  quality  portion  of 
the  model  simulates  temperature  and  eight  water  quality 
constituents  including  dissolved  oxygen  and  phytoplankton.  The 
model  will  determine  the  water  quality  needed  from  all  reservoir 
releases  to  meet  specified  downstream  water  quality  objectives 
and  will  determine  the  gate  openings  in  each  reservoir  that  will 
yield  the  appropriate  reservoir  release  water  quality.  Should  it 
be  necessary,  flows  will  be  altered  to  ensure  that  downstream 
water  quality  objectives  are  met.  The  model  selects  the  "best" 
solution  for  the  system-wide  reservoir  operation. 


CHAPTER  III 


RESERVOIR  SYSTEM  DESCRIPTION 

The  Sacramento  Valley  reservoir  system  consists  of  four 
major  reservoirs  as  shown  in  Figure  4.  Shasta  and  Keswick  are 
tandem  reservoirs  in  parallel  with  Oroville  and  Folsom 
Reservoirs.  Numerous  tributaries  and  irrigation  diversions  are 
involved . 

Shasta  and  Keswick  Reservoirs  are  located  on  the  Sacramento 
River  in  northern  California  in  the  northern  end  of  the 
Sacramento  Valley  about  240  river  miles  north  of  Sacramento. 
Keswick  is  a  reregulation  reservoir  designed  to  even-out  the 
daily  hydropower  releases  from  Shasta.  Below  Shasta  and  above 
Keswick,  inter-basin  water  transfers  enter  the  Sacramento  River 
through  Spring  Creek.  Along  the  Sacramento  River,  Cow  Creek  and 
Cottonwood  Creek  are  major  inflowing  tributaries  and  the 
Anderson-Cottonwood ,  Tehama-Colusa,  Corning  and  Glenn-Colusa 
Irrigation  District  Canals  are  major  irrigation  diversions. 

Oroville  Reservoir  is  located  on  the  Feather  River  in  the 
Sierra  foothills  about  95  river  miles  north  of  Sacramento.  Major 
tributaries  entering  the  Feather  River  include  the  Yuba  and  Bear 
Rivers.  Major  diversions  are  located  immediately  below  Oroville 
Dam  from  the  Thermalito  Afterbay.  The  Feather  River  flows  into 
the  Sacramento  River  near  Verona. 

Folsom  Reservoir  is  located  on  the  American  River  in  the 
Sierra  foothills  about  30  miles  east  of  Sacramento.  The  American 
River  below  Folsom  Reservoir  is  leveed  with  no  major  tributaries 
entering  before  its  confluence  with  the  Sacramento  River  at 
Sacramento  (I  Street) . 


The  Sacramento  River  continues  to  flow  south  towards  the  San 
Francisco  Bay  via  the  Sacramento  and  San  Joaquin  Delta.  This 
study's  lower  boundary  is  located  near  Hood  about  20  miles  south 
of  Sacramento  and  just  upstream  of  the  Delta's  maze  of 
interconnected  waterways . 
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CHAPTER  IV 


APPLICATION  PROCEDURE  . 

i 

The  application  of  the  HEC-5  to  the  Sacramento  Valley 
reservoir  system,  described  in  Chapter  III,  or  to  any  other  \ 

system  begins  with  data  collection,  assembly,  and  manipulation.  » 

Each  of  these  tasks  in  addition  to  model  execution  is  discussed  J 

below.  Interpretation  of  results  is  discussed  in  Chapter  V.  t 

i 

Data  Collection 

The  HEC-5  model  data  requirements  are  similar  to  those  of 
most  comprehensive  water  quality  models.  The  data  to  be 
collected  are  categorized  into  three  types;  time  independent, 
required  time  dependent  and  optional  time  dependent. 

The  time  independent  data  include:  physical  description  of  j 

the  reservoir  (i.e.,  elevation  vs.  volume,  surface  area  and  dis¬ 
charge  capacity;  and  vertical  reservoir  segmentation),  physical 
description  of  the  river  (i.e.,  river  mi?e  vs.  cross  section  and 
channel  discharge  capacity;  and  river  reach  segmentation) ,  con¬ 
trol  point  desired  and  required  flows,  model  coefficients  (i.e., 
flow  routing;  reservoir  diffusion;  physical,  chemical  and  biolo¬ 
gical  reactions  rates)  and  initial  conditions  for  the  start  of 
the  simulation.  The  input  data  for  the  Sacramento  Valley  reser¬ 
voir  system  is  shown  in  Appendix  A.  The  time  independent  data 
for  the  water  quantity  input  are  on  the  RL,  RO,  RS,  RQ,  RA,  RE, 

R2 ,  CP,  ID  and  RT  records.  For  the  water  quality  input,  the  time 
independent  data  are  on  LI,  L2,  LR,  L3,  L6,  L7,  L8,  PL,  L9,  Cl, 

C2,  C5 ,  C6 ,  C7 ,  SA,  DK,  CR,  SI,  S2 ,  SR,  S3,  S4  and  KR  records. 

The  required  time  dependent  data  include:  evaporation,  wind 
speed,  cloud  cover,  air  temperature,  dew  point  temperature,  flow 
diversions,  inflow  quantity  and  quality  for  all  reservoir  and 
river  tributaries,  discharge  quantity  from  reservoirs,  and 
control  point  target  water  quality  conditions.  The  time 
dependent  data  are  on  the  R3,  IN,  QA  and  QD  records  for  the  water 
quantity  input.  The  input  for  the  water  quality  portion  are  on 
the  EZ,  ET,  CT,  12  and  14  records. 

The  optional  time  dependent  data  include:  reservoir 
storages;  river  flows  at  other  than  control  points;  and  reservoir 
water  quality  profiles  and  river  time  series  plots.  These  data 
are  used  as  checks  (using  auxiliary  graphics  programs)  on  the 
model  output  in  contrast  to  the  previously  mentioned  data  which 
are  required  to  make  the  model  work.  This  data  is  referred  to  as 
calibration  data. 

Data  Assembly 

Sources  for  the  data  categorized  above  are  numerous.  In 
general,  they  include  all  water-related  agencies  at  the  federal, 
state,  local  and  private  levels.  To  name  a  few,  the  following 
should  all  be  considered  when  searching  for  data: 


t 

f 


Corps  of  Engineers  (COE) 

Bureau  of  Reclamation  (BUR) 

Geological  Survey  (USGS) 

Environmental  Protection  Agency  (EPA) 
state  departments  of  natural  resources 
state  environmental  protection  agencies 
state  colleges  and  universities 
city  and  county  public  works  offices 
utility,  water  and  flood  control  districts 


Meteorological  data  (dry  bulb  temperature,  dew  point 
temperature,  cloud  cover  and  wind  speed)  are  readily  available 
from  the  U.S.  Weather  Service,  local  airports  and  universities. 
The  primary  data  source  for  Corps  of  Engineers  offices  is  the 
U.S.  Air  Force  Environmental  Technical  Applications  Center  (OL-A 
USAFETAC)  and  for  other  offices  is  the  NOAA's  National  Weather 
Service  (NWS)  office;  both  in  Asheville,  North  Carolina.  They 
should  be  contacted  early  in  the  data  assembly  task,  and  an  order 
should  be  placed  for  the  required  meterological  data.  Normally 
it  will  take  30  to  60  days  to  receive  the  data  after  the  request 
is  made.  This  data  service  is  free  to  the  Corps  of  Engineers  and 
is  documented  in  Department  of  Army  Pamphlet  115-1.  It  is 
recommended  that  data  format  CD  144  be  requested  from  OL-A 
USAFETAC  at  (FTS)  672-0218  or  NWS  at  (704)  259-0682. 


Data  Manipulation 


For  simplification,  the  meteorological  data  should  be 
reformatted  with  the  "Weather"  computer  program  [HEC  1986]  and 
then  input  to  the  "Heat  Exchange"  computer  program  [U.S.  Army 
1977],  available  from  HEC,  which  provides  output  in  the  form 
necessary  for  the  HEC-5  program  (i.e.,  ET  record  images).  If 
HEC-5  computation  intervals  shorter  than  daily  are  going  to  be 
used,  the  use  of  the  "Heat  Exchange"  program  should  be 
substituted  with  a  similar  utility  program. 


The  cross-sectional  data  may  be  already  available  in  the 
proper  format  if  previous  flood  plain  or  flood  damage  analysis 
were  studied  using  the  HEC-2  program.  Because  HEC-2  is  a  widely 
used  water  surface  profile  program,  the  author  has  found  the 
cross-sectional  data  already  available  in  the  proper  format  on 
many  past  water  quality  studies.  For  simplification,  the  HEC-2 
type  data  should  be  input  to  the  "GEDA"  computer  program  [HEC 
1981],  available  from  the  HEC,  which  provides  output  in  the  form 
necessary  for  the  HEC-5Q  program  (i.e.,  S3  card  images). 


The  meteorological  and  cross-sectional  data  account  for 
about  60%  of  the  water  quality  portion  of  the  data  to  be  prepared 
on  a  typical  study. 


Tributary  inflows,  diversions  and  reservoir  discharges  may 
be  readily  available  from  WATSTORE  and  STORET  data  systems. 
WATSTORE  is  managed  by  the  USGS  and  contains  streamflow  data. 
STORET  is  managed  by  the  EPA  and  contains  water  quality  data. 
These  computer  data  systems  can  often  provide  the  necessary 
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inflow  and  reservoir  discharge  quantity  and  quality  data. 
Modification  of  output  format  to  the  required  HEC-5  input  format 
for  the  IN,  QA,  QD  and  14  cards  is  trivial.  Past  experience  has 
shown  that  the  flow  data  from  WATSTORE  is  often  not  available  at 
the  same  location  as  water  quality  data  from  STORET.  Therefore 
judgement  must  be  used  to  determine  its  usability  by  direct 
transfer  of  location  or  appropriate  adjustment.  The  two  types  of 
data  must  be  compatible  for  the  computer  program  to  calculate  the 
appropriate  constituent  load  time  series.  This  type  of  data  is 
about  20%  of  the  water  quality  portion  of  the  input  data. 

The  remaining  20%  of  the  typical  data  set  is  located  by 
searching  through  COE  reservoir  regulation  manuals,  searching 
through  COE  files  and  making  numerous  judgement  decisions 
regarding  a  thorough  knowledge  of  the  study  objectives,  the 
drainage  area  being  studied  and  the  concepts  employed  by  the 
HEC-5  model. 


The  model  can  be  used  in  several  different  ways.  The 
cheapest,  fastest  execution,  also  requiring  the  least  data 
preparation,  is  a  steady  state  analysis.  The  input,  calculation 
and  output  intervals  are  monthly.  These  may  be  used  for 
screening  monthly  data  for  multiple  years  to  find  a  critical 
period  (poor  water  quality  condition)  for  more  detailed  analysis. 
This  method  will  only  find  the  most  critical  period  of  at  least 
monthly  duration  and  not  necessarily  critical  periods  of  shorter 
duration. 

The  model  can  be  used  to  study  temperature  only  or 
temperature  and  any  combination  of  up  to  three  conservative 
parameters  (i.e.,  chlorides,  total  dissolved  solids  -  TDS, 
alkalinity,  specific  conductance) ,  up  to  three  nonconservative 
parameters  (i.e.,  coliform,  carbonaceous  biochemical  oxygen 
demand  -  CBOD,  ammonia)  and  dissolved  oxygen  (DO) .  If  the 
phytoplankton  option  is  requested,  the  parameters  in  the  model 
must  include  at  least  temperature,  TDS,  nitrate,  phosphate, 
phytoplankton,  CBOD,  ammonia,  and  DO.  The  more  parameters 
simulated,  the  more  computation  and  data  preparation  time 
involved.  The  phytoplankton  option  is  particularly  time 
consuming  since  all  the  parameter  calculations  possible  are  being 
performed  with  this  option. 

The  model  can  be  used  for  existing  and/or  proposed 
reservoirs.  If  an  existing  condition  is  being  simulated,  usually 
the  objective  is  to  reproduce  historical  events  through  model 
calibration.  The  calibration  option  can  take  about  50%  less 
computer  time  because  the  time-consuming  linear  programming 
algorithms  are  not  used. 

Calibration  should  begin  by  reproducing  observed  temperature 
profiles  in  the  reservoirs  and  stream  channels.  The  reservoir 
diffusion  coefficients,  A1  and  GSWH,  provide  the  best  initial 
adjustments.  The  second  step  in  reservoir  calibration  involves 
adjustments  to  the  three  factors  affecting  light  penetration: 


secchi  disk;  solar  radiation  absorbed  near  surface,  XQPCT;  and 
depth  associated  with  XQPCT  amount  of  solar  radiation.  The  above 
five  variables  are  interrelated  and  will  need  to  be  adjusted 
simultaneously.  Although  many  other  variables  affect  the  thermal 
and  water  quality  constituent  calibration  of  reservoirs,  it  is 
not  recommended  to  adjust  them  unless  you  have  more  data  than 
normal  for  reproducing  observed  profiles.  Reservoir  calibration 
decisions  may  include  the  user's  choice  of  the  input  weather 
station. 

The  calib  ation  of  the  river  profiles  may  also  involve  the 
users  choice  of  weather  stations  to  index  the  meteorology 
affecting  the  river.  Unless  more  than  normal  data  are  available, 
it  is  not  recommended  to  adjust  any  model  variables. 

Once  the  model  has  been  calibrated,  the  objective  may  be  to 
modify  an  existing  reservoir  operation  pattern  or  to  evaluate  the 
impact  of  proposed  new  reservoirs  or  channel  modifications.  This 
analysis  requires  the  use  of  the  linear  programming  algorithm  and 

the  increase  of  computer  tine  is  significant. 

The  simulation  mode  discussed  above  can  be  used  either  to 
evaluate  the  best  water  quality  that  can  be  provided  throughout 
the  system  for  given  reservoir  discharges  (obtained  either 
external  to  the  simulation  or  determined  by  the  HEC-5  quantity 
part  of  the  model)  or  to  evaluate  the  best  water  quality 
operation  without  preconceived  discharge  quantities.  The  former 
operation  is  referred  to  as  a  balanced  pool  operation  and  the 
latter  as  a  flow  augmentation  operation. 

The  balanced  pool  operation  is  the  standard  HEC-5  analysis 
for  flow.  When  using  the  balanced  pool  operation,  the  water 
quality  portion  of  the  program  simply  evaluates  the  best  vertical 
level  for  withdrawal  at  each  reservoir  (assuming  multiple  level 
intakes  are  available)  to  meet  all  downstream  water  quality 
targets  for  the  given  reservoir  discharge. 

The  flow  augmentation  operation  allows  the  model  to  relax 
the  balanced  pool  concept  and  to  decide  how  much  flow  should  come 
from  which  reservoir  and  at  which  vertical  level  in  order  to  meet 
downstream  water  quality  targets.  Sometimes  downstream  water 
quality  improvements  require  significantly  increased  discharge 
rates  to  obtain  only  small  improvements  in  water  quality.  This 
flow  augmentation  operation  is  the  most-costly  execution  and  not 
always  a  practical  alternative  for  real  world  regulation  given 
power,  water  conservation  and  flood  control  storage 
considerations . 

For  this  demonstration  application,  the  input  data  shown  in 
Appendix  A  was  executed  using  the  calibration  option. 

Application  of  this  option  allows  the  user  to  define  the  exact 
level  of  the  intake  structure  operated  and  the  exact  quantity  of 
discharge  from  each  dam.  This  is  the  normal  method  of  model 
application  when  calibrating  the  model  to  observe^  uisuoricai 
data.  This  application  was  executed  for  water  temperature, 
specific  conductance,  alkalinity,  CBOD,  ammonia,  and  dissolved 
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oxygen  because  this  was  the  data  available.  Actually,  all  these 
types  of  data  are  very  limited  in  availability  on  most 
tributaries  but  some  of  them  (e.g.,  temperature,  dissolved  oxygen 
and  specific  conductance)  are  more  readily  available  as 
in-channel  data.  The  in-channel  data  are  the  optional  time 
dependent  data  used  for  calibration. 

Selected  portions  of  the  computer  output  are  shown  in 
Appendix  B.  Graphical  displays  of  the  results  are  shown  in 
Appendix  C  and  are  discussed  in  Chapter  V. 
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7. 


CHAPTER  V 
SIMULATION  RESULTS 

The  Sacramento  Valley  reservoir  system  operation  described 
in  Chapter  III  was  simulated  using  HEC-5  and  produced  results 
which  were  compared  to  water  quantity  and  quality  data  in  the 
four  reservoirs  and  at  all  downstream  control  points.  The  data 
for  comparison  purposes  consisted  of  discharge  rates  at  most 
control  points  as  well  as  water  temperature  at  many  of  the  same 
locations.  Other  water  quality  parameters  are  less  available  but 
are  compared  where  they  are  available. 

The  model  simulation  for  the  Sacramento  Valley  system  used 
temperature,  specific  conductance  (sometimes  called  electrical 
conductivity) ,  alkalinity,  carbonaceous  biochemical  oxygen  demand 
(CBOD) ,  ammonia  (NH3)  and  dissolved  oxygen  (DO) .  These  specific 
parameters  were  chosen  based  on  the  availability  of  at  least 
limited  data  except  for  CBOD  which  was  estimated  and  adjusted  by 
calibration  to  reproduce  DO. 

The  Sacramento  Valley  reservoir  system  results  are  shown  in 
Appendix  B  in  an  abbreviated  form.  The  computer  output  begins 
with  an  echo  of  the  input  data.  These  output  should  be  examined 
carefully  to  insure  that  the  program  is  getting  the  data  that  the 
user  expects  the  program  to  execute.  This  step  of  interpretation 
of  results  is  important  and  well  worth  the  required  time. 

The  remainder  of  the  computer  output  are  the  day-by-day 
results  of  water  quality  profiles  in  the  reservoirs  and  along  the 
stream  network.  The  output  is  quite  voluminous  and  only  selected 
portions  have  been  included.  Because  it  is  so  voluminous,  it  is 
also  difficult  to  interpret  in  tabular  form. 

The  graphical  display  of  these  results  is  included  as 
Appendix  C  for  the  reservoirs  and  at  selected  locations  along  the 
stream  network.  While  this  graphical  capability  is  not  a  part  of 
HEC-5,  an  option  exists  for  the  model  to  write  the  results  to  an 
HEC  data  storage  system  called  DSS  [HEC  1985].  An  HEc  graphics 
program  called  DSPLAY  (described  in  the  DSS  manual)  can  read  the 
results  from  DSS  and  produce  plots  as  shown.  These  plots 
satisfactorily  demonstrate  the  capability  of  HEC-5  to  reasonably 
reproduce  observed  reservoir  and  stream  profiles  on  large 
systems . 

The  legends  at  the  bottom  of  the  reservoir  water  quality 
plots,  pages  C1-C16,  define  simulated  and  observed  data  for 
various  dates.  Shasta,  Oroville  and  Folsom  Reservoirs  have 
sufficient  observed  temperature  data  to  be  useful  for  calibration 
purposes.  Sufficient  observed  data  for  the  other  parameters  were 
not  available  except  in  Folsom  Reservoir  as  shown. 
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Considering  the  model  limitation  of  having  only  one  weather 
station  for  the  entire  system  (subsequently  modified  to  use 
parameter  statements),  it  is  the  author’s  opinion  that  the 
reproduction  is  quite  good.  Perhaps  some  further  refinement 
could  be  achieved  with  additional  trials,  but  the  acceptability 
of  the  model  can  be  demonstrated  with  these  results.  The 
Oroville  Reservoir  observed  data  anomaly  in  June  (page  C6)  cannot 
be  explained  with  the  available  input,  particularly  when  the  much 
cooler  observed  hypolimnion  temperatures  of  May  and  July  are 
noted. 

The  legend  at  the  bottom  of  the  time  series  stream  plots 
defines  the  various  observed  and  simulated  water  quality 
parameters  for  the  study  period.  As  with  the  reservoir  plots, 
only  those  locations  which  have  sufficient  observed  data  to  be 
used  for  calibration  purposes  are  shown.  Unlike  the  simulated 
data,  the  observed  data  points  are  often  more  than  one-day  apart. 
Extreme  caution  should  be  applied  to  any  interpolation  between 
observed  data  points. 

In  general,  the  calibration  of  the  model  is  quite  good  along 
the  Sacramento  River  for  all  the  observed  parameters  downstream 
to  Hamilton  City.  Butte  City  and  Colusa  temperatures  (pages  C24 
and  C26  respectively)  show  that  significant  warming  of  this  reach 
of  the  Sacramento  River  takes  place  at  least  during  the  Spring 
(April  and  May  1979) .  This  temperature  consideration,  in 
addition  to  the  lack  of  sufficient  simulated  quantity  of  flow  at 
Butte  City  and  Colusa  (pages  C23  and  C25  respectively)  suggests 
that  the  undefined  return  flows  on  the  Sacramento  River  between 
Hamilton  City  to  Knights  Landing  are  sufficiently  large  to  cause 
significant  errors.  Other  parameter  reproductions  are  also  poor; 
apparently  due  to  undefined  return  flows. 

The  Feather  River  below  Oroville  and  the  American  River 
below  Folsom  lack  sufficient  water  quality  data  to  provide 
adequate  information  for  calibration  purposes.  The  reproduction 
of  observed  flow  is  shown  in  pages  C29  and  C30. 

Careful  interpretation  and  evaluation  of  the  Sacramento 
River  results  lead  the  author  to  encourage  the  continued 
application  of  this  model  to  help  develop  understanding  of  the 
workings  and  operation  of  any  stream  system. 


CHAPTER  VI 


SUMMARY 

HEC-5  can  be  used  for  comprehensive  water  quality  studies 
involving  complex  river  network  and  reservoir  systems.  The 
program  is  used  to  compute  the  best  operation  for  a  reservoir 
system  and  determine  either  the  water  quality  condition  resulting 
from  the  best  water  quantity  operation  (balanced  pool  method)  or 
the  best  water  quality  operation  without  maintaining  balanced 
reservoir  conservation  pools  for  the  system  (flow  alteration 
method) . 

The  HEC-5  model  uses  a  linear  optimization  scheme  to 
determine  the  target  water  quality  at  the  dam  to  best  meet  all 

user  weighted  downstream  targets.  Then  a  non-linear  scheme  is 

used  with  user  weights  to  determine  how  the  intake  structure  will 
be  operated. 

The  Sacramento  River  results  and  results  of  the  two 
additional  applications  referenced  in  Chapter  II  have  each 
provided  model  improvements  and  added  confidence  in  the  model 
validation.  While  it  is  true  that  large  comprehensive  data  sets 
have  not  been  used  to  validate  the  model,  the  author  believes 

that  the  variety  of  years  modeled  and  the  variety  of  locations 

studied  have  provided  sufficient  experience  to  the  Corps  of 
Engineers  to  warrant  the  continued  use  and  application  of  the 
HEC-5  water  quality  model. 

It  is  the  conclusion  of  this  study  that  HEC-5  is  a  viable 
tool  for  evaluating  reservoir  systems  operation  for  water  quality 
analysis. 
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HEC-5  Input  for  Sacramento  River  System 
Water  Quality  Modeling 
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